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Abstract 

In fermion mixing phenomenology, the matrix of moduli squared, P = {\U\ 2 ), 
is well-known to carry essentially the same information as the complex mixing 
matrix U itself, but with the advantage of being phase-convention independent. 
The matrix K (analogous to the Jarlskog CP-invariant J) formed from the real 
parts of the mixing matrix "plaquette" products is similarly invariant. In this 
paper, the P and K matrices are shown to be entirely equivalent, both being 
directly related (in the leptonic case) to the observable, locally L/E'-averaged 
transition probabilities in neutrino oscillations. We study an (over-)complete set 
of flavour-symmetric Jarlskog-invariant functions of mass-matrix commutators, 
rewriting them simply as moment-transforms of such (real) invariant matrices. 
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1. Introduction: The P matrix 

The complex mixing matrix (CKM matrix [1] or MNS matrix [2] respectively) occupies 
pride of place at the center of quark and lepton mixing phenomenology. The mixing 
matrix connects mass-eigenstates of different charges, for example, charged-leptons to 
neutrinos, through an array of complex mixing amplitudes: 
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where our choice to display the mixing matrix for leptons here, rather than quarks, 
anticipates our emphasis on neutrino oscillations below. Of course, the freedom to 
rephase rows and columns of the mixing matrix means that the mixing-matrix elements 
(in this case the U\ v where I — e, //, r and v — 1, 2, 3) are phase-convention dependent. 

While no one would deny the importance of amplitudes and the mixing-matrix 
concept, it has long been appreciated [3] that observables (including the magnitudes 
of CP-violating asymmetries 1 see below) are always ultimately expressible in terms 
of moduli-squared of mixing elements, whereby we often find it convenient to work 
instead with the corresponding matrix of moduli-squared, P = {\U\ 2 ) — UoU* (where 
the "o" denotes the simple entrywise (Schur) product of two matrices): 





\U e l\ 


2 


\Ue2\ 


2 


\Ue3\ 




p = (\U\ 2 ) = UoU* = ^ 




2 




2 




i) 




u Tl \ 


2 


u T2 \ 


2 


u T3 \ 





The P-matrix represents rather the probabilistic content of the various states (again, in 
this case, the charged-letpton mass-eigenstates versus the neutrino mass-eigenstates) 
and has the advantage of being independent of any choice of phase convention. 

The P-matrix for leptons (Eq. 2) features naturally in neutrino oscillation phe- 
nomenology. We have, eg. that P P T (where T denotes the matrix transpose) repre- 
sents directly [4] the flavour-to-flavour transition/oscillation-probability matrix V(L/E) 

<V>oo := lim <V(L/E)> = PP T (3) 

L/E-^oo 

for neutrino oscillations in the "asymptotic" domain (oo), ie. when L/E » (Am-,)" 1 
for all neutrino mass-squared differences Amf^, i,j — 1 — 3, with local L/P-averaging 

1 Of course the information on the sign of the relevant Jarlskogian J (for quarks or leptons) is lost in 
taking the moduli, with a consequent overall sign ambiguity affecting all associated CP-asymmctries 
(for quarks or leptons respectively). In principal, however, this single qbit of information (per sector) 
could be carried in an ad hoc fashion by giving a sign to the P-matrix itself, which will anyway be 
taken as positive in any other context. We stress that we expect the considerations of this paper to 
be valid for both Dirac and Majorana neutrinos, provided we restrict attention to flavour oscillation 
phenomena, neglecting neutrinoless double beta-decay observables. 
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(L the propagation length and E the neutrino energy). The corresponding result for 
the "intermediate" domain (N), when (AmJ 2 ) _1 >> L/E » (Am| 3 ) _1 , (Am^)" 1 , 
again averaging over unresolved oscillations within the domain, is also expressible 
entirely in terms of the P-matrix (see below). The transition/oscillation probabil- 
ity matrix V(L/E) gives disappearance probabilities on the diagonal and appearance 
probabilities off the diagonal, and (given the sign of J) its full L/E- dependent form, ie. 
with no L/E'-averaging, is also expressible in terms of the P-matrix, see Eq. 8 below. 

We begin by noting that any two rows or columns of the P-matrix are sufficient to 
calculate the magnitude of the Jarlskog CP- violation parameter [5] using, respectively: 



4J 2 
4J 2 



(Pll P V 1 + Pl2 Pl<2 + PlS Pl'3 f - 2 (f* Py x + Pi P? 2 + Pi P? 3 ) (4) 

\2 o/D2 d2 , D2 r>2 , r>2 r>2 



{PevPev' + PfivP^v' + PtvPtv' ) ~ 2 [P ev P ev > + F ' -F ' / + P TU P TV ,) (5) 



where / 7^ /' and v 7^ u' . Symmetrising over, eg. rows, and exploiting the Schur product 
and the 3x3 identity matrix / together with its binary complement / 
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we obtain a relatively simple matrix expression for J 2 : 
1 



J 



24 



Tr[(/o (PP 7 ))(/o {PP 1 )) -21(1 o ((PoP)(PoP) i ))] 



(7) 



which is quartic in P and obviously P <-> P T symmetric. 

The full time-dependent flavour-to-flavour vacuum transition/oscillation probabil- 
ity matrix V(L/E) is ('dimensionally') quadratic in P and may now be written: 



V(L/E) = Io (P 



+ 7o(P 
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1 (S12 + S23 + S31) (8) 
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where the first term gives the disappearance probabilities and the second and third 
terms the CP-even and CP-odd parts of the appearance probabilities (the latter con- 
tributing proportional to the 'epsilon matrix' e [6]). The time dependence in Eq. 8 
enters only through c^- := cos Am|L/2£ and Sij := sin Am^L/2E. The asymptotic 
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limit is obtained setting <c i j> OQ — <s i j> 00 — — 1 — 3) when Eq. 8 reduces 
to PP T as expected. The corresponding probabilities in the intermediate domain are 
obtained setting <c 12 >n = 1, <si2>« = <C23>m = <s 2 3>m = <C3i>n = <S3i>n = 0. 

We may remark that average appearance probaility measurements in the interme- 
diate domain determine directly the z/ 3 column of the P-matrix, via <V e ^>^= 2P e3 P fl3 , 
<V^ t >m= 2P fl3 P r3 etc. whereby P e3 = [<P eM >M<P er >N / <V^ T >^ etc. (a zero in 
the P-matrix here leads to indeterminacies, in which case we must use disappearance 
probabilities also: <P^> N = 1-2P m3 + 2P 2 3 P^ = l/2±[l/4-(l- <? w > M )/2]5). 
Disappearance probabilities in the asymptotic domain then determine the remaing in- 
dependent column: <V ee > 00 = P^ + P e 2 2 + P e 2 3 , <V^>oo = P*i + Pfe + P% etc. 
whereby P el(2) = (1 - P e3 )/2 ± [(1 - P e3 ) 2 /4 + P e3 (l - P e3 ) - (1- <P ee >oo)/2]^ etc. 
(sign ambiguities correlate with the interchange of entries between P-matrix columns). 
Note that a complete set of average vacuum probability measurements, in the inter- 
mediate and asymptotic domains together, in general over-determines the P-matrix 
and so fixes the magnitude of J (via Eq. 7), despite the fact that the CP- violating 
third term in Eq. 8 makes no contribution to averaged probabilities in either domain. 

Clearly we must fore-go here any detailed discussion of experimental data, measure- 
ment errors, matter effects etc. Indeed it will suffice, for the purposes of this paper, 
simply to take as given, in the first instance, the succinct summary of the current 
oscillation data provided by the tri-bimaximal ansatz [7] [8]: 



'2/3 1/VS \ /2/3 1/3 

U = | _i /v /6 1/V3 1/V2 \ =► P = V6 1/3 1/2 ] (9) 
-1/VQ 1/V3 -I/V2J \l/6 1/3 1/2, 

where we display the mixing matrix and the corresponding P-matrix side-by-side. 
From, eg. Eq. 4, applied to, eg. the first two rows of the P-matrix Eq. 9, we recover: 
j2 = ((2/3)(l/6) + (l/3)(l/3)) 2 - 2((2/3) 2 (l/6) 2 + (l/3) 2 (l/3) 2 ) = Q (1Q) 

Comparing the LHS and RHS of Eq. 9, the P-matrix definitely provides the more 
visually clear summary, without the distracting arbitrary (unobservable) phases/signs. 

It turns out that we shall also find it useful below to refer, occasionally and for 
illustrative purposes, to the older trimaximal ansatz [9] (now ruled out by the data): 

(l/y/3 1/a/3 1/-n/3\ /1/3 1/3 1/3 \ 

U= lu/VS 1/V3 u/V3\ =>• P= 1/3 1/3 1/3 (11) 

\lu/V3 1/V3 u/VS) \l/3 1/3 1/3/ 

{uj = exp(i2-7r/3) and cu = exp(—i2n/3) represent the complex cube-roots of unity). 
From Eq. 4/Eq. 5, applied to any two rows/columns of the P-matrix Eq. 11, we have: 

j2 = ((1/3) 2 + (1/3) 2 + (1/3) 2 ) 2 - 2((l/3) 4 + (1/3) 4 + (1/3) 4 ) = J_ 

4 108 1 1 
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corresponding, as expected, to maximal CP-violation (J = ±l/6\/3 [10]). 

The extreme hierarchical form of the CKM matrix, ie. the smallness of the off- 
diagonal elements, makes the P-matrix in the quark case much less useful visually, eg. 
as a summry of data, than it is in the leptonic case. Since in this paper we shall be 
emphasising the P-matrix and in particular some useful matrices computable from it 
(see Sections 2-4) we shall henceforth focus essentially exclusively on the leptons. 

2. Plaquettes: The K matrix 

A mixing matrix "plaquette" [11] is a 4-subset of complex mixing matrix elements, 
obtained by deleting one row and one column of the mixing matrix. Associating a 
given plaquette with the element common to the deleted row and column, there is 
clearly one plaquette 'complementary to' each mixing element. The corresponding 
plaquette product may be defined: IL^ := U\-\ u -iU*_i v+1 Ui + i u +iU* + i v _\ (generation 
indices to be interpreted mod 3, meaning that when / ="r" we have I + 1 ="e" etc.). 
The matrix of plaquette products (for leptons) may then be written: 

(u t3 u* t2 u» 2 u; 3 u Tl u; z u^u;, u^u^u;^ 
n = u e3 u* e2 u T2 u; 3 u el u: 3 u r3 u* Tl u e2 u* el u Tl u; 2 . (13) 

\U, 3 U; 2 U e2 U* e3 U^UaUZ U^U*JJ el U^J 
Plaquette products are phase-convention independent complex numbers, and, as is 
well-known [12], all have the same 2 imaginary part, equal to and defining the Jarlskog 
CP-violating parameter J [5]. The real parts define the .fT-matrix (H^ := —Ki u + iJ): 

(K el K e2 K e3 \ (J J J\ 

n = — K^i K, 2 7^3 + i\J J J (14) 
\K Tl K t2 K t3 J \J J J) 

(where for our convenience we have introduced a minus sign in the definition with 
respect to some previous work, K Xv := — Kfy" [14\) . The if-matrix [13] [14] may be 
viewed as the natural CP-conserving analogue of the Jarlskog variable J. 

The i^-matrix is an important oscillation observable. The magnitnude of the Jarl- 
skog CP violation parameter may be obtained directly [13] from the i^-matrix by 
summing products in pairs of K-matrix elements within any row or column: 

K n K l2 + K l2 K l3 + K l3 K n = J 2 (15) 
K ev K^ + K^K TU + K TU K eu = J 2 . (16) 

Symmetrising over rows or columns we obtain respectively: 

J 2 = ^Tr [ IKK T ] = ^Tr [ IK T K ] (17) 



2 Very often [12] alternating signs are found to enter here (so that the imaginary parts in Eq. 14 
are given by ±J). However, with our cyclic definition (Eq. 13), there are no such alterning signs. 
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which is clearly very succinct and obviously K <-> K T symmetric. 

The full time-dependent flavour-to-flavour vacuum transition/oscillation probabil- 
ity matrix is ('dimensionally') linear in K and may be written: 



V(L/E) = Io(I-2IK 



+1+ o (2I_K 



1 1 - C 23 1 - C 2 3 1 - C 23 
1 - C31 1 - C31 1 - C31 
V 1 - C12 1 - C12 1 - C12 




C23 
C31 
Ci2 



1 - C 23 N 
1 - C31 
1 - Ci2 / 



) 



\S\2 + «23 + S 3 1 



where we exploit the (cyclic) flavour raising and lowering operators: 
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(19) 



(J = I + + J_ and e = i+ — J_). Note that, in Eq. 18, the second and third terms, 
constituting the CP-even part of the appearance probabilities, are simply transposes 
one of the other. The asymptotic limit may be written: 



<V> OQ = I + (/+ - /) o (2I.KD) + (/_ - /) o (2I+KD) 



(20) 



where D is the so-called 'democratic matrix' (D = I + I), and appearance probabilities 
are now given by the if-matrix row-sums: <V e ^>oo = 2(if T i + K r2 + K r3 ) etc. In the 
intermediate domain, appearance probabilities are given by: <P eM >M = 2(if Tl + K t2 ) 
etc. so that the z/ 3 column is determined directly: K t3 = (<P eM >oo — <V e ^>M)/2 
etc. and J 2 may be calculated from Eq. 16. The remaining elements are given by: 
K Tl{2 ) = <V e ^>M /4±[<V efl >i /16+ <V e ^>M K t3 /2-AJ 2 }^, where the sign ambiguity 
correlates with the interchange of the v x and v 2 columns, as before. 

The i^-matrx is readily computable [13] [14] from the P-matrix (again, indices to 
be interpreted mod 3) using: 



K lv = {P iv - {Pi 



Pi. 



l-l u-l r l+l u+l 



p 



Pi. 



l-l u+l r l+l u-l 



0)/2, 



(21) 



ie. just take the corresponding element of the P-matrix, subtract the permanent of its 
associated "P plaquette" (defined analogously to a mixing matrix plaquette above), 
and divide by an overall factor of two (the permanent of a matrix is defined similarly 
to its determinant, but with the alternating signs replaced by all positive signs). 
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Given instead the i\"-matrix we can also always obtain the P-matrix, in effect 
inverting Eq. 21. Specifically, to obtain the P matrix element P^, we focus on the cor- 
responding K matrix element K iv and especially on its complementary "i\"-plaquette" 
(again defined analogously to a mixing matrix plaquette). We sum the products of 
adjacent elements of the if-plaquette giving A\ v and add twice its permanent, M tl/ : 

A lv = (K l+1 „ +1 + Ki^i v ^i)(Ki-i „ +1 + K l+1 (22) 
Mi„ = Ki +1 u+ iKi-x v -\ + Ki-i u+ iKi + i v -\ (23) 

{Ai v is the product of the sums of diagonally opposing elements of the -fT-plaquette 
while M lv is the sum of their products). Finally we divide by the sum of all K-elements 
not in the if-plaquette (excluding K\ v itself) and take the (positive) square root: 3 

Pi „ = ^Av + 2M lu )/{K l+1 u + KiZ u + K x u+1 + Kx Zij. (24) 

To illustrate the above procedures we can take the tri-bimaximal ansatz as a con- 
venient example, and display the P matrix and the K matrix side- by-side.: 

2/3 1/3 \ /1/6 1/12 -1/18\ 

1/6 1/3 1/2 K = 1/9 . (25) 

1/6 1/3 1/2/ V V 9 / 

In computing K from P we have, eg. that the K e \ element is given by P e \ less its 
complementary P-permanent, P^PtI + P^PtZ, all divided by two: 

K el = (2/3 - (1/3 x 1/2 + 1/2 x l/3))/2 = (2/3 - l/3)/2 = 1/6. (26) 

In computing P from K we have, eg. that for P ei the .^-permanent is zero, and there 
is only one non-zero adjacent-pair product (P^3 x P T 3), whereby: 



Pel 



1/9 x 1/9 



1/12 - 1/18 \ 1/36 



1/9 x 1/9 




(27) 



with only two terms {K e2 = 1/12, K e3 = —1/18) non-zero in the denominator sum. 
The Jarlskog parameter computed from the K- matrix Eq. 25, eg. from the first row: 

is seen to vanish, as expected for tribimaximal mixing. 



3 Since the argument of this square root must always be positive for a valid K- matrix, we could 
equally (c.f. footnote 1) carry the sign of the CP- violation adopting an ad-hoc sign convention for K. 
(For non-zero mixing, the sum of all if-matrix elements is always positive in our present convention.) 
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Similarly (and for further illustration) we have for the trimaximal ansatz: 

/1/3 1/3 l/3\ / 1/18 1/18 1/18\ 

P= 1/3 1/3 1/3 K= 1/18 1/18 1/18 (29) 

V 1/3 1/3 1/3/ \1/18 1/18 1/18/ 

In computing K from P we have that, eg. K e \ is given by: 

K el = (1/3 - (1/3 x 1/3 + 1/3 x l/3))/2 = (1/3 - 2/9)/2 = 1/18 (30) 
and in computing P from K we have, eg. for P el : 



Pel = 



if (18 x 18) 



4/(18 x 18) + 2 x 2/(18 x 18) 
\J 4/18 ~ \J 4/l8 




(31) 



where clearly all elements of either matrix are equivalent in the trimaximal case. From 
any row or column of the i^-matrix, summing products of pairs of elements, the 
Jarlskog invariant is given by J 2 = 3/18 2 = 1/(36 x 3), ie. J = ±1/(6^3) as expected. 

So if we have the full if-matrix we can calculate the full P-matrix and vice-versa. 
The i^-matrix and the P-matrix are thus entirely equivalent, both carrying complete 
information about the mixing, including the CP violation (excepting only its sign). 
In Appendix A we show that, just as the P-matrix is determined by any P-plaquette, 
the .fT-matrix is similarly determined by any if-plaquette (up to a twofold ambiguity). 

3. Quadratic Commutator Invariants: The Q-matrix 

The P, K and U matrices above all have rows and columns (respectively) labelled 
by the charged-lepton and neutrino mass-eigenstates. For some applications, however, 
flavour-symmetric mixing variables, making no reference whatever to any particular 
basis, could be considered more appropriate (e.g. as input to the proposed "extremisa- 
tion" programme [15]). The prototype for such variables is undoubtedly the Jarlskog 
CP- violating invariant itself [5], written here for the leptonic case in terms of the com- 
mutator C := — i[L, N] of the charged-lepton (L) and neutrino (N) mass matrices: 

TrC 3 /3= -2 Detdiag(Ai) Det diag(A„) J := -2L A N A J (32) 
Aj = (m M - m T , m T - m e , m e - m M ), = (m 2 - m 3 , m 3 - mi, mi - m 2 ) (33) 

with m e , m M , m T the charged-lepton masses and m 1; m 2 , m 3 the neutrino masses. 
In this section and the final section respectively, we shall examine, systematically, 
generalisations of the Jarlskog invariant, to a range of quadratic and cubic commutator 
traces, which turn out to be intimately related to the K and P matrices (Sections 1-2). 
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We begin by considering quadratic commutator invariants [16] defined (C„ 
-i[L m ,N n ] [16]) as follows [15]: 



Q = 



l / TrC 2 u TrC7 n C7 12 TrC 2 12 \ 

~2 ^ r( ^ 11< ^ 21 TrCnC22 TrCi2C22 , 
\ TrC^ TTC21C22 Tr C| 2 / 



(34) 



each trace clearly flavour-symmetric and vanishing in the case of zero mixing. Row 
and column labels now correspond to successive powers of charged-lepton and neutrino 
mass matrices respectively Note that powers higher than L 2 or N 2 need not be consid- 
ered, since, eg. L 3 can be always re-expressed in terms of L 2 , L and L° = I, by virtue 
of the characteristic equation. Furthermore Tr [L 2 , N] [L, N 2 } = Tr [L, N] [L 2 , N 2 }, so 
that there are indeed just nine such invariants. 

One may then show, eg. by use of the flavour projection operator technique [17], 
that the Q nm (n,m = 1, 2, 3) are just suitably-defined double-moments of the if-matrix: 



Q nm = A/ diag(A / ) (diagE,) K (diagS,)™" 1 diag(A„) A, 



(35) 



with mass-sum vectors £/, E„ (cf. the mass-difference vectors A/, A u above) given by: 

£j = ("fy + m r, ~m T + m e , m e + m M ) E„ = (m 2 + m 3 , m 3 + mi, m x + m 2 ). (36) 

Eq. 35 (taken together with the definition Eq. 34) is the manifest analogue of Eq, 32 for 
the CP-conserving case. While the Q-matrix is then simply a linear transform of the 
.fT-matrix, with respect to appropriate lepton and neutrino Vandermonde matrices: 



Q 



( 1 

(m M + m T ) 
\ (m M + m T f 

K (diagA,) 2 



{m T + m e ) (m e + m M ) | (diagA;) 2 

(m T + m e ) 2 (m e + m M ) : 

1 1 1 

(m 2 + m 3 ) (m 3 + mi) (mi + m 2 ) 

(m 2 + m 3 ) 2 (m 3 + mi) 2 (mi + m 2 ) 2 



(37) 



the i^-matrix is clearly simply the corresponding inverse transform of the Q-matrix: 
(diagA; 



K = 



1 I (m e + m^)(m T + m e ) -(2m e + m M + m r ) 1 
(m^ -\- m T )(m e + m^) — (2m M + m r + m e ) 1 
^ A \ (m T + m e )(m M + m T ) — (2m r + m e + m M ) 1 

(mi + m 2 )(m 3 + mi) -(2mi+m 2 + m 3 ) l\ T ^ ; „„ A 
Q I (m 2 + m 3 )(mi + m 2 ) -(2m 2 + m 3 + mi) 1 



(diagA^)" 



(38) 



(m 3 + m 1 )(m 2 + m 3 ) -(2m 3 +m 1 +m 2 ) 1 
so that (for known masses) the Q-matrix and the K-matrix are entirely equivalent. 
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Of course i^-elements are not all independent; products of pairs of elements in any 
row or column sum to J 2 [13] [14]. Corresponding relations result for the Q-matrix: 

P?-i^H-i(Qii*3i3 — Qn) + ^-1^+1(^-1 + Si + i)(2Qi2<522 — Q11Q23 — Q13Q21) 

+ + £; +1 )(2<522Q32 ~~ Q21Q33 — Q23Q31) + (Q31Q33 — Q32) 

+S / _ 1 S i+1 (Q 11 Q 3 3 + Q13Q31 - 2Q 12 Q32 + 4Q21Q23 - 4Q 2 ! 2 ) 

+A 2 (Q 21 Q 23 - Q 2 2 )]/(A Z L A iV A ) 2 = J 2 (39) 

(Q11Q31 - Q21) + + ^Wl) (2Q21Q22 — Q11Q32 ~ Q31Q12) 

+ + S iy+ i)(2Q22<523 — Q12Q33 — Q32Q13) + (Q13Q33 — Q23) 

+Tj p _iTj v+ i{QhQ^ + Q13Q31 ~~ 2Q21Q23 + 4Q12Q32 — 4Q22) 

+A 2 (Q 12 Q 32 - gL)]/(A,L A iV A ) 2 = J 2 (40) 

for I — e, /i, t, and z/ = 1, 2, 3, respectively (where for eg. / = e, = S M = m T + m e , 
= S T = m e + m^, A/ = A e = — m T etc.). Equating such expressions among 
themselves yields relations between Q-matrix elements, and in particular flavour- 
symmetric relations may be obtained, eg. by summing over all / and all v respectively: 

N\[L 2 PA {Q\iQiz — Q\ 2 ) + 2L p^L p 3 (2Qi 2 Q 22 — Q11Q23 — Q13Q21) 

+ 2Lp 3 Lp 2 (2Q 22 Q 32 — Q21Q33 ~ Q23Q31) + -^P 2 (Q3lQ33 ~~ Q32) 

+L P6 (QnQ33 + Q13Q31 - 2Q12Q32 + 4Q 2 iQ 23 - 4Q 2 ! 2 ) + 3L A Q 21 Q 23 ] 
= L 2 A [N 2 P4 (Q n Q3i ~ Q21) + 2N P4 N P3 (2Q 21 Q22 - Q11Q32 - Q31Q12) 
+2N P3 Np2(2Q 2 2Q23 ~ Q12Q33 ~ Q32Q13) + N 2 P2 (Q 13 Q 33 - Q 2 23 ) 
+N P6 (Q n Q 33 + Q 13 Q 31 - 2Q 21 Q 23 + 4Q 12 Q 32 - 4Q 2 2 ) + 3iV 2 Q 12 Q 32 ] (41) 

(basically the relation Tr IKK T = Tr IK T K). Further flavour-symmetric relations 
could clearly be obtaned, eg. by equating the sums over products in pairs etc., but 
such higher-power relations would not seem to have much practical value here. Clearly 
only five relations among Q-matrix elements can be functionally independent in total. 

The supplementary polynomials used in Eq. 41 above are defined in terms of the 
traces of powers of mass matrices: L x := Tr L, L 2 : = Tr L 2 , Aq := Tr N etc. [15]. 
For, eg. the charged-leptons, we have: 

L E := {L\ - L 3 )/3 = (m e + m M ) (m^ + m T ) (m r + m e ) (42) 

L P2 : = (3L 2 - Ll)/2 = m 2 e + m 2 + m 2 - m e m^ - m^m T - m T m e (43) 
L P3 := (3L 3 -L 2 L 1 )/2 = m3 + mJ + m^ 

-(m 2 m M + m 2 m e )/2 - (m 2 m T + m 2 m M )/2 - (m 2 m e + m 2 m T )/2 (44) 
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L P4 := L 2 P2 + 2(LiL 3 - L\) = + + ra* — m^m 2 ^ - m^m 2 . - m 2 T m 2 e (45) 

L 2 A := A(Lp A L P2 - L 2 P3 )/3 = (m e - m M ) 2 (m M - m T ) 2 (m T - m e ) 2 (46) 
L P6 := L 2 P3 - L 2 A / A = (ml + m 3 ^ + ml - m^m^ - m 2 ^m T - mlm e ) 

x (m e + + m r — m e m^ — m p jn T — m T m e ) (47) 

with analogous relations for the neutrino polynomials in terms of N ± , N 2 , N 3 . 

4. Cubic Commutator Invariants: P and R Matrix Moments 

Finally, we consider all cubic commutator trace invariants, completely generalising 
the original Jarlskog CP-invariant (as for the mass- matrices themselves, powers of 
commutators higher than cubic can aways be reduced via the characteristic equation). 

Taking into account the cyclic property of the trace we may classify the cubic 
commutator invariants, according to the number of repeated commutators entering. 
For example, there are four such invariants with all three comutators identical: 



-iTrQ\ = 6iL A N A J (48) 

-iTrCf 2 = 6iNxL A N A J (49) 

-iTr C\ x = 6iLxL A N A J (50) 

-iTrC 2 3 2 = 6iL^Nj]L A N A J (51) 



all clearly proportional to J, and carrying no information on mixing angles beside the 
value of J itself, assuming the masses are known. 

Then one should consider the twelve cubic commutator invariants with two com- 
muatators identical and one different: 



-*Tr C 2 U C 12 


= AiN\L A N A J 


(52) 


-iTr C 2 U C 21 


= 4iL 1 L A N A J 


(53) 


-iTr C 2 2 C n 


= 2iNp 2 L A N A J 


(54) 


-iTi C 2 2l C u 


= 2iLp 2 L A N A J 


(55) 


—iTr C 2 2 C 22 


= 4iN^LiL A N A J 


(56) 


-iTr C\ X C 22 


= AiL-£,NiL A N A J 


(57) 


—iTr C\ 2 C\ 2 


= 2iN^L P2 L A N A J 


(58) 


-*Tr C 2 22 C 2l 


= 2iL^N P2 L A N A J 


(59) 


-iTi C 2 1X C 22 


= 2i(T 11 + L 1 N 1 )L A N A J 


(60) 


-iTi C 2 12 C 21 


= 2i(— T 12 + L 1 N 2 )L A N A J 


(61) 


-tTr C 2 2l C l2 


= 2i(-T 21 + L\N 1 )L A N A J 


(62) 


-iTr C 2 2 2 Cn 


= 2%(T 22 + L\Nl - (L\ + L 2 ){Nl + N 2 )/A)L A N A J 


(63) 
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The first eight of these (Eqs. 52-59) have no dependence on mixing angles except 
through J as before, and together with Eq. 48 above would in fact be more than suffi- 
cient to fix all the masses in terms of cubic commutators alone, if that was desirable. 
The last four traces however (Eqs. 60-63) having no triply repeated index in either 
sector, are seen to involve double- moments (T mn [12]) of the P- matrix: 

T mn = Tr L m N n = mf(diag m,)™ -1 P (diag m v ) n - x m v (64) 
mj = (m e ,m ll ,m T ) ml = (m 1 , m 2 , m 3 ) (65) 

with the moment- weightings here defined in terms of the simple mass- vectors mi, m„. 
Clearly T 00 = 3, T 01 = N 1: T 02 = N 2 etc. by unitarity. The lowest four non-trivial 
moments T n , T 12 , T 21 ., T 22 appearing in Eqs. 60-63 above, are clearly sufficient to 
determine the P-matrix itself, assuming the masses are known. 

Finally, we consider the cubic commutator invariants with all three commutators 
different which are seen to be complex (the R m ,n are defined just below Eq. 74): 

-iTr C U C 12 C 21 = R ll L A N A + i{-T 11 + ^L 1 N 1 )L A N A J (66) 

-iTr C 12 C 22 C U = R 12 L A N A + i(Ti 2 + Li(2iVf - N 2 ))L A N A J (67) 

-iTr C 21 C U C 22 = R 21 L A N A + i(T 21 + (2L 2 1 -L 2 )N 1 )L A N A J (68) 
— iTr C 22 C 2 \C\ 2 = R 22 L A N A + 

i{-T 22 + {2>L\Nl - L\N 2 - NlL 2 + N 2 L 2 )/2)L A N A J. (69) 

Readily constructed from the P and K matrices, the R- matrix is real and given by 4 : 

Riv = (Pi-iu-iKi-iu — Pi-iu+\Ki u+ i + Pi + i u+ \Ki + i u — Pi + i u ^\Ki u ^i). (73) 

For example, the i?-matrix for tribimaximal mixing is given by: 



'2/3 1/3 

1/6 1/3 1/2 | => /,' 
1/6 1/3 1/2 



/ 0\ 

-1/18 1/18 . (74) 

V i/i8 -i/i8 o; 

4 The Riv are also expressible in terms of differences of "hexaplaquettes" (which are mixing in- 
variants comprising products of six mixing matrix elements, choosing two elements (one complex- 
conjugated) from each row and column, as noted by Jarlskog and Kleppe [18]). We have: 

&f v ~ (U l ^ 1 „U?_ lu+1 Ui iy+1 Ui„_ 1 U l+liy - 1 U? +lu ) 

-{Ui-iv-iUUvUiv+xUtv-M+ivUt+^y*/* (70) 
:= -((^-i v +in, v+ i + fl + i v _in? v _j-(fl_i v _in,_i v + ^ + i v+ in; +lw )**/*) (7i) 

:= -Ri„-i{{Pi-i v+ i- Pi+iv-i)T{Pi-iv-i- Pi+iv+i))J (72) 

where the sign ambiguity correlates with the absence or presence of a relative complex conjuga- 
tion in the difference, indicated in Eq. 70 by the overall * * /* appearing on the second product. 
Hexaplaquettes feature for example in calculations of the electric dipole moment of the neutron [19]. 
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In Eqs. 66-69 (RHS), the real parts are proportional to the P-matrix moments, which 
are written simply R mn (R mn := mf (diag m;) m_1 P (diag m I/ ) n ~ 1 m I/ , cf. Eq. 64). The 
associated imaginary parts are proportional to J, but depend linearly also on the T mn . 
Note that the P-matrix for trimaximal mixing is simply the null matrix, identical in 
fact to the P-matrix for no-mixing at all, whereby it is already clear that we cannot 
unambiguously recover the P and K matrices starting from the P-matrix in general. 

Considering cubic and quadratic traces (Section 3) together, we clearly have here 
a functionally (over-)complete set of flavour-symmetric Jarlskog-invariant mixing vari- 
ables, simply related to the P and K matrices (Sections 1-2). Assuming that flavour 
symmetry and Jarlskog invariance are always to be respected [20], we may anticipate 
that our present analysis will prove pertinent to formulating future theories of flavour. 
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Appendix A 

For the P-matrix we know that we need only the four elements of any plaquette 
to obtain all the elements, since any row or column of the P-matrix sums to unity. 
Analogously, for the A-matrix we can calculate the full matrix (at least up to a twofold 
ambiguity, see below) in terms of the four elements of any A-plaquette. We need the 
sum {Si v ) of the elements of the A-plaquette, the product of all four elements (Fi v ), 
and the sum of the products taken in threes (Wi„) and also the determinant (Di u ): 

Siu = Ki+i u+i + Ki-i „+i + Ki + i „_! + Ki^i „_! (75) 

Fiu = Ki+i u+iKi-i u+iKi+i u-iKi-i u-i (76) 

W lv = F lv ( 1/AVi „+i + „+i + u-i + 1/AVx „_! ) (77) 

Diu = Ki+i v+iKi-i v-i — Ki-i u+iKi+i v-i (78) 

where the formula for W\ v has been written here in terms of F\ v only for brevity (ie. 
in the case of any A-plaquette elements being zero, it is to be understood that those 
denominators are to be cancelled with the corresponding factor in the numerator before 
proceeding with the substitution.) The square of the Jarlskog invariant J is given by: 

/ 2 - ( Wlv + Dl )( 1 + f 1 9 i Dl m\ 



13 



The ^-elements not in the K-plaquette (with the exception of the K x v element itself) 
can then normally be obtained using, eg. : 

Ki+i v = (J 2 - K l+ i u+1 Ki +1 „_i)/ (Ki+i v+1 + K l+ i (80) 

with analogous formulae for K\_\ v , K il/+1 etc. The K lv element itself is given by: 

Kx „ = -(J 4 - (A lu + P lv )J 2 + F lu )/(S lu J 2 - W lv ). (81) 

As an example of these procedures we shall first consider the i^-matrix for tri- 
maximal mixing Eq. 29 (since the tribimaximal case requires a little more more care). 
Suppose that only the K r3 plaquette is known, ie. we have K ei = 1/18, K e2 = 1/18, 
K^2 = 1/18, Kpi = 1/18, but no a priori knowledge of K e3 , K^ 3 , K Tl , K t2 or K r3 . We 
wish to reconstruct the full trimaximal K matrix: 



1/18 1/18 *' 
K = \ 1/18 1/18 * | => K 



^ 1/18 1/18 1/18\ 
1/18 1/18 1/18 (82) 
1/18 1/18 1/18/ 



From the K r3 plaquette (Eq. 82 LHS) we have: 

^ = 4/18 F r3 = l/18 4 W^ t3 = 4/18 3 D t3 = (83) 
We may then calculate J 2 from Eq. 79 as follows: 



J 2 = (1/(4 x 18) + 0)(l/2 ± - 4/18) - (84) 
'' x(-±-) = — 1— or — 1— (85) 



4x18 v 2 6 y 4 x 27 8 x 27 

Taking J 2 = 1/(4x27) corresponding to maximal CP violation (J = J max = ±1/6a/3), 
ie. corresponding to trimaximal mixing, we have, eg. 

1/(4 x 27) - 1/18 x 1/18 1/(2x81) 
Ke3 = 1/18 + 1/18 = 2/18 = 1/18 (86) 



and 



1/(4 x 27) 2 - (4/18 2 + 2/18 2 ))/(4 x 27) + 1/18 4 
r3 4/18 x 1/(4 x 27) -4/18 3 [ ' 

(9 - 18 + l)/(2 4 x 3 8 ) 1 



(3-l)/(2x3 6 ) 18 

which leads to Eq. 82. Taking instead the alternative solution J 2 = 1/(8 x 27) (with 
intermediate CP violation J = J max /V2 = ±l/6\/6) leads to the K matrix (RHS): 

/5/12 5/12 l/6\ /1/18 1/18 1/72 \ 

P= 5/12 5/12 1/6 K= \ 1/18 1/18 1/72 (89) 

V 1/6 1/6 2/3/ \l/72 1/72 23/144/ 
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corresponding to the uninteresting P matrix also displayed (LHS). Thus the K matrix 
is determined from a K plaquette here, only up to a two-fold ambiguity. 

In case of zeros in the denominators above, we must use less succinct but slightly 
more general formulae (as in the definition of Wi v ). The element K i+ i „ then becomes: 



Ki+i u = 



Bi+i v , 1 . / 1 



Ki+i u+\Ki + i v -i(Ki-i v+ i + Ki-i „_i) + K l+1 u+ iKf_ x u _ l + K l+1 u ^iKf_ x v 



+1 



02 



(90) 



where B i+1 v is the product of three sums of adjacent pairs of elements as follows: 

Bi+i v = (Ki + i v+ i + K X -\ v+1 )(Ki + i v _i + K X -\ v+ i + Ki-i „_i) (91) 

ie. excluding the pair in line with the element required. The Ki u element is given by: 

2S lv G lv + H lv (l±y/l-4S lv ) 



K lv 



(92) 



Sl{2S lv -{l±y/l-AS lv )) 
where G\ v and H\ v are defined by: 

G lv := 8S lv F lv + S lv P lv (S? v -2M lv ) + 3DlS lv -Dl-W lv S lv (l + S lv ) (93) 
H lv := W lv S lv -Sf v (A lv + M lv ) + (2S lv -l)(AF lv -M^). (94) 

Tri-bimaximal mixing is a case in point. Suppose we wish to reconstruct the tri- 
bimaximal K matrix, starting only from its K r3 plaquette: 



/1/6 1/12 *\ 
K= * 

in which case we have: 



K 



/1/6 1/12 -1/18 \ 

1/9 
V 1/9 / 



A t3 = 1/72 ^3 = 1/4 
From Eq. 79 the CP violation is zero (J 2 

K e3 ~- 

We have then also: 



(95) 



(96) 



1/6 x 1/12 



M t3 = F t3 = W t3 = D t3 = 0. 

: 0) and, eg. the K e3 element is given by: 
4 1 



1/4 



6 x 12 



18' 



G r3 = 

so that the K r3 element is given by: 



H r3 = -1/4 3 x 1/72. 



K, 



e3 



-1/4 3 x 1/72 _ 4 x 2 _ 1 

~ 9 



l/4 4 x(-l/2) 72 
with no ambiguities arising (clearly no ambiguitities arise whenever Si„ = 1/4). 



(97) 
(98) 

(99) 
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